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INTRODUCTION 

During  the  past  decade,  the  Navy  has  considera¬ 
bly  improved  its  capabilities  in  aircraft/weapon 
integration.  In  1989  it  took  more  than  400  hours 
of  wind  tunnel  testing,  which  cost  1,5000,000, 
and  20  flights1  to  clear  the  JSOW  from  the  F-18 
to  Mach  0.95.  This  year  the  MK-83  JDAM  was 
cleared  after  only  60  hours  of  wind  tunnel  testing 
and  five  flights  to  the  full  F-18  aircraft  envelope 
of  Mach  1.3.  This  reduction  occurred  because 
the  Navy  not  only  learned  to  test  smarter,  but  also 
developed  an  integrated  approach2  to  Modeling 
&  Simulation  (M&S),  wind  tunnel  and  flight 
testing  which  allowed  lessons  learned  on  previ¬ 
ous  programs  to  be  applied  to  new  ones. 
However,  the  present  approach  still  requires  a 
fairly  large  commitment  of  time  and  financial  re¬ 
sources  to  accomplish  the  mission. 

The  present  approach  has  optimized  the  use 
of  available  resources;  any  further  gains  will  have 
to  come  not  by  improving  existing  techniques, 
but  by  bringing  new  resources  into  the  process. 

Several  years  ago  the  Office  of  the  Secretary 
of  Defense  (OSD),  under  the  Central  Test  and 
Evaluation  Investment  Program  (CTEIP)  funded 
a  tri-service  research  project  termed  Applied 
Computational  Fluid  Dynamics  (ACFD)  for  store 
separation.  This  project  is  meant  to  provide 
analysis  tools  that  effectively  use  Computational 
Fluid  Dynamics  (CFD)  for  store  certification 
analysis.  ACFD  will  provide  the  needed  tools 
that  will  reduce  DOD  dependence  on  wind  tunnel 
and  flight-testing. 


ACFD  is  not  intended  to  replace  the  wind 
tunnel  in  the  near  future;  rather  it  will  be  used  to 
determine  the  critical  regions  of  the  flight  enve¬ 
lope  to  help  structure  the  wind  tunnel  test,  and  to 
explain  any  wind  tunnel  anomalies  and  help 
structure  the  flight  test  program.  The  objective 
of  the  program  is  to  provide  upgraded  analysis 
tools  that  will  support  store  certification  require¬ 
ments  at  less  cost  and  in  less  time. 

NOMENCLATURE 

BL:  Aircraft  Buttline,  positive  outboard,  in. 

Cp  Rolling  moment  coefficient,  rt  wing  down 
Cm:  Pitching  moment  coefficient,  positive  up 
Cjq:  Normal  Force  coefficient,  up 
Cn:  Yawing  moment  coefficient,  nose  right 
Cy:  Side  force  coefficient,  right 
FS:  Aircraft  Fuselage  Station,  positive  aft,  in. 
l.e.  F-18  Wing  leading  edge  flap 
t.e.  F-18  Wing  trailing  edge  flap 
M:  Mach  number 

P:  Store  roll  rate,  positive  rt  wing  down 

Q:  Store  pitch  rate,  positive  nose  up 

R:  Store  yaw  rate,  positive  nose  right 

PHI:  Store  roll  angle,  positive  rt  wing  down, 
deg. 

PSI:  Store  yaw  angle,  positive  nose  right,  deg. 
THE:  Store  pitch  angle,  positive  nose  up,  deg. 
WL:  Aircraft  Waterline,  positive  up,  in. 

Z:  Store  C.G.  location,  positive  down,  ft. 

a:  Angle  of  attack,  deg. 

e:  Upwash  angle,  positive  up,  deg. 
ct:  Sidewash  angle,  positive  outboard,  deg. 

Note:  all  wind  tunnel  data  shown  are  right  wing, 
flight  test  left  (negative  PSI,  PHI,  Y) 
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ACFD  CHALLANGES 

Over  the  past  several  years  there  have  been  sev¬ 
eral  organized  efforts  to  validate,  demonstrate 
and  accelerate  the  insertion  of  CFD  methods  into 
the  store  certification  process  for  external  stores 
carriage  and  release.  Several  significant  efforts 
have  been  documented  in  AIAA  conference  pro¬ 
ceedings.  The  first  of  these  was  the 
Wing/Pylon/Finned- Store  which  occurred  in 
Hilton  Head,  SC  in  the  summer  of  1992.  An  ex¬ 
tensive  set  of  wind  tunnel  store  carriage  and 
separation  data  for  CFD  code  validation  were 
made  available  for  a  generic  wing  and  store  ge¬ 
ometry1 2 3 *.  Although  Euler4,5  and  thin  layer  Navier 
Stokes6  (TNS)  solutions  were  in  good  agreement 
with  these  test  data,  solution  times  on  the  order 

H 

of  5  days  on  the  Cray  YMP  made  such  tools  im¬ 
practical  for  everyday  use.  Madson  later 
demonstrated8  that  the  TranAir  full-potential 
code  could  give  results  of  similar  quality  in  a 
fraction  of  the  time  required  for  the  higher  order 
codes,  as  may  be  seen  in  Figure  1  and  Figure  2. 


WING  UPPER  SURFACE 

M  =  0.95  BL  =  5.3 


LOCATION  X/C 

FIGURE  1 


WING  LOWER  SURFACE 

M  =  0.95  BL  =  5.3 


ACFD  CHALLENGE  I 

The  first  ACFD  sponsored  conference  was  the  F- 
16/Generic  Finned  Store9'14  which  occurred  in 
New  Orleans  in  the  summer  of  1996  (ACFD 
Challenge  I).  At  the  end  of  the  meeting,  the 
ACFD  tri- service  technical  leads  evaluated  the 
CFD  tools  that  were  used  to  predict  the  F-16  Ge¬ 
neric  store  carriage  loads.  The  evaluation 
concentrated  on  the  following  characteristics: 

1.  Time  required  to  obtain  a  solution,  which  in¬ 
cluded  the  time  required  to  convert  the  geometry 
into  a  form  compatible  with  the  code’s  pre¬ 
processor,  the  time  required  to  obtain  a  surface 
and  volume  grid,  and  the  time  required  to  ensure 
that  the  solution  obtained  was  properly  con¬ 
verged. 

2.  The  accuracy  of  the  solution. 

3.  Code  efficiency  or  time  required  to  obtain  a 

meaningful  minimum  number  of  independent 


LOCATION  X/C 

FIGURE  2 


solutions;  one  solution,  at  one  Mach  number  and 
angle  of  attack  is  useless  for  determining  a 
store’s  trajectory.  This  time  is  determined  both 
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by  the  difficulty  in  generating  a  new  geometry 
and  the  running  time  for  each  case. 

4.  Potential  for  production  use:  the  overall  goal 
of  the  ACFD  project. 

The  following  CFD  codes  were  evaluated 
for  the  project: 

SPLITFLOW9  by  Lockheed. 

TranAir10  by  NAWCAD 

USM3D12  by  NAWCAD 

OVERFLOW13  by  NASA  Johnson 

STORESIM14  by  Wright  Labs 
Since  STORESIM  never  obtained  a  valid  solu¬ 
tion,  it  was  dropped  from  the  study. 

I.  The  time  required  to  obtain  a  solution,  ac¬ 
cording  to  the  reports  submitted  at  the  end  of  the 
project,  were  hard  to  compare  since  it  could  be 
influenced  to  a  large  extent  by  the  past  experi¬ 
ence  of  the  code’s  user.  The  NAWCAD  efforts 
were  purposefully  conducted  by  individuals  with 
no  prior  experience  with  the  particular  codes 
used,  while  the  other  efforts  were  conducted  by 
individuals  with  varying  degrees  of  prior  experi¬ 
ence  with  the  tools.  However,  based  to  the 
information  available,  TranAir  took  about  two 
weeks  to  obtain  a  solution,  SPLITFLOW  and 
USM3D  took  about  a  month  for  a  solution,  while 
OVERFLOW  took  about  three  months  for  one 
solution.  On  that  basis  the  TranAir  code  was 
judged  the  best. 

n.  The  accuracy  of  the  solution  was  also  hard  to 
determine,  since  SPLITFLOW  had  24  independ¬ 
ent  solutions,  some  of  which  were  excellent 
while  others  were  terrible,  while  OVERFLOW 
had  only  one.  However,  based  on  the  one  com¬ 
mon  solution  for  the  two  codes,  it  appears  that 
the  OVERFLOW  solution  was  slightly  better  in 
terms  of  the  pressure  comparisons.  The  USM3D 
and  TranAir  solutions  were  comparable  to  each 
other,  but  not  as  good  as  the  OVERFLOW  com¬ 
parison. 

III.  Since  SPLITFLOW  provided  24  solutions  in 
the  four  months  of  the  project,  TranAir  and 
USM3D  provided  8  independent  solutions,  and 
OVERFLOW  had  only  one  solution,  the 


SPLITFLOW  code  was  clearly  superior  in  this 
respect. 

IV.  For  transition  potential  SPLITFLOW  ap¬ 
pears  to  be  clearly  superior  to  all  the  others. 
TranAir,  which  gives  answers  of  equivalent  accu¬ 
racy  with  considerably  less  set  up  time  seemed 
slightly  better  than  USM3D.  OVERFLOW,  as 
the  Navy  had  previously  determined,  is  too  cum¬ 
bersome  to  provide  meaningful  results  in  the  time 
required  for  store  separation  projects. 

Many  important  lessons  were  learned;  how¬ 
ever,  the  experimental  test  case  did  not  include 
flight  test  data  ("real"  store  trajectories).  Because 
of  this  limitation,  store  certification  engineers 
continued  to  express  skepticism  towards  the  ac¬ 
curacy  of  CFD  methods.  Also,  the  CFD 
community  raised  concerns  about  the  credibility 
of  portions  of  the  wind  tunnel  test  data,  criticiz¬ 
ing  scale,  model  support  interference,  and  wall 
effects.  Therefore,  there  was  a  desire  within  the 
ACFD15  program  to  reconcile  these  issues  by 
conducting  additional  analysis  by  using  a  data  set 
that  included  both  wind  tunnel  and  flight  test 
data. 

ACFD  CHALLENGE  II 
Selection  of  Test  Case 

Large  sets  of  wind  tunnel  and  flight  test  data  ex¬ 
isted  for  the  F/A-18C  JDAM  configuration. 
Figure  3.  During  the  flight  test  phase,  both  pho- 
togrametrics  and  telemetry  were  used  to  track  the 
position  of  the  store  during  releases.  Out  of  these 
tests,  two  release  conditions  were  selected  for 
this  CFD  Challenge.  The  basis  for  these  two 
cases  included  the  following  considerations:  1) 
matching  aircraft  and  store  geometry  in  both 
wind  tunnel  and  flight  tests,  2)  correlation  be¬ 
tween  wind  tunnel  data  and  flight  test  data,  3) 
possession  of  both  high  transonic  and  low  super¬ 
sonic  cases  with  interesting  miss  distance  time 
histories,  4)  ability  to  publicly  release  the  wind 
tunnel  and  flight  test  data  to  an  international 
audience. 
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FIGURE  3  F/A-18C/JDAM  Test  Case 


TEST  CASE  PARAMETERS 

The  test  cases  selected  were  M  =  0.962  at  6,382 
ft.  (flight  13)  and  M  =  1.05  at  10,832  ft.  (flight 
14).  Both  cases  were  for  the  aircraft  in  a  45- 
degree  dive. 

For  these  two  test  cases,  the  configuration 
geometry  for  the  wind  tunnel  and  flight  test  is 
shown  above.  The  JDAM  is  mounted  on  the 
outboard  pylon,  with  the  330-gallon  fuel  tank  on 
the  inboard  pylon.  The  SUU-65  BRU-32A/A 
ejector  rack  provided  a  nominal  peak  force  of 
7,000#  for  both  fwd  and  aft  cartridges.  As  may 
be  seen  in  Figure  4,  an  Euler  prediction  of  the 
Mach  number  distribution  at  M  =0.962  indicates 
strong  transonic  interference  effects. 

FLIGHT  TEST  RESULTS 

Both  Captive  Trajectory  System  (CTS)  grid  data, 
and  store  aerodynamic  force  and  moment  data 
measured  on  the  wing  pylon  were  available  for 
this  aircraft  configuration.  These  data  were  input 
into  a  six-degree-of-freedom  trajectory  code  be¬ 
fore  the  flight  tests  were  performed.  Parametric 
variations  on  flight  conditions  and  store  aerody¬ 
namic  forces  were  performed  to  ensure  that  the 
flight  test  could  be  safely  accomplished.  After 
the  flight  tests  were  completed,  the  trajectory 
simulations  were  again  performed,  with  the  ac¬ 
tual  flight  conditions  used  to  try  to  match  the 
flight  test  results  These  predicted  trajectories 
were  used  a  the  metric  for  the  CFD  Challenge  II, 


FIGURE  4  JDAM  MK-84 


and  were  not  given  to  the  participants  until  just 
prior  to  the  meeting 

TEST  FLIGHT  #13 

Flight  test  #13  was  conducted  on  July  10,  1996. 
The  store  was  released  in  a  43  degree  dive  at 
6,382  ft.  at  M  =  0.962.  The  telemetry  and  photo- 
grametric  data  were  not  in  good  agreement  with 
each  other  for  the  vertical  displacement.  Since 
inertial  effects  (store  mass  and  ejector  force) 
largely  drive  the  vertical  displacement,  the  rela¬ 
tive  Z  displacement  is  usually  the  easiest  to 
predict.  The  discrepancy  in  Z  was  attributed  to 

JDAM  FLIGHT  13 
M  =  0.962  6382  FT  43  DIVE 


DEG 


FIGURE  5 
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the  effects  of  aircraft  motion  caused  by  store  re¬ 
lease. 

As  may  be  seen  in  Figure  5,  the  predicted 
pitch  and  yaw  attitudes  at  M  =  0.962  were  in  ex¬ 
cellent  agreement  with  the  flight  test  results.  The 
roll  attitude  was  not  well  predicted. 

However,  roll  attitude,  which  is  the  hardest 
to  predict,  fortunately  has  a  minimal  impact  on 
the  trajectory.  The  photogrametric  results  are  not 
shown,  since  they  were  considered  to  be  less  ac¬ 
curate  than  the  telemetry  data. 

TEST  FLIGHT  #14 

Flight  test  #14  was  conducted  on  August  29, 
1996.  The  store  was  released  in  a  44  degree  dive 
at  10,832  ft.  at  M  =  1.055.  The  telemetry  and 
photogrametric  results  for  the  displacement  again 
showed  a  large  discrepancy  in  Z. 

The  predictions  using  the  wind  tunnel  test 
data  were  again  in  excellent  agreement  with  the 
flight  test  data,  Figure  6. 


MISS  DISTANCES 

The  most  important  parameter  for  safe  store 
separation  is  the  ability  to  accurately  predict  the 
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FIGURE  7  JDAM  MISS  DISTANCE 

JDAM  FLIGHT  14 
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FIGURE  8  JDAM  MISS  DISTANCE 
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store  miss  distance,  which  is  the  smallest  dis¬ 
tance  between  any  part  of  the  store  and  aircraft 
during  the  early  part  of  the  trajectory.  The  tra¬ 
jectory  simulation  used  the  force  measured  for 
ejection  from  the  F-18  centerline  to  compensate 
for  aircraft  motion. 

As  may  be  seen  in  Figures  7  and  8,  the  miss 
distance  predictions  were  in  excellent  agreement 
with  the  test  data.  The  disagreement  between  the 
photogrametric  and  telemetry  predicted  miss 
distance  is  attributed  to  the  fact  that  the  telemetry 
could  not  take  aircraft  motion  into  account. 

ACFD  CHALLENGE  II  PAPERS 

General 

Each  participant  was  requested  to  include  in  his 
or  her  paper: 

1)  a  description  of  the  CFD  and  trajectory  inte¬ 
gration  methods  used  to  produce  the 
estimates  of  the  trajectory; 

2)  a  description  of  the  methods  and  resources 
required  to  produce  the  computational  grid; 

3)  estimates  of  carriage  loads,  the  position  and 
attitude  of  the  store  throughout  the  computed 
trajectories  and  an  estimate  of  the  miss  dis¬ 
tance  versus  time;  and 

4)  metrics  of  the  CFD  process  used,  including 
convergence  rate,  man-hours  and  time  re¬ 
quired  for  grid  generation,  computer 
resources  used  and  an  estimate  of  the  exper¬ 
tise  of  personnel  required  to  replicate  the 
results. 

Eight  papers  by  Cenko16  ,  Hall17,  Tomaro18, 
Woodson19,  Welterlen20,  McGroy21,  Fairlie22,  and 
Benmeddour23  were  submitted  for  ACFD  Chal¬ 
lenge  II.  The  meeting  was  held  at  the  AIAA 
Annual  meeting  in  Reno  Nevada  on  January  12th, 
1999.  Due  to  the  interest  in  the  Challenge  the 
timing  of  the  session  and  the  venue  were  changed 
to  enable  seating  for  around  200  people;  despite 
this,  the  room  was  filled  to  capacity  with  over  50 
people  having  to  stand  in  the  back  for  four  hours. 


The  first  paper16  described  the  wind  tunnel 
and  flight  test  results,  while  the  other  seven  de¬ 
scribed  the  application  of  seven  different  CFD 
codes  to  the  problem.  Two  of  the  papers17,21 
were  not  ready  in  time  to  be  included  in  the 
meeting  proceedings,  but  all  eight  papers  were 
either  presented  at  the  meeting,  or  the  results 
were  provided  at  a  later  date. 

ACFD  Challenge  II  Overview  paper 

Cenko  outlined  the  background  to  the  Challenge 
and  the  sources  of  the  data  used.  The  wind  tun¬ 
nel  data  consisted  of  both  CTS  grid  and  carriage 
force  and  moment  data  measured  on  the  wing 
pylon,  conducted  in  the  CALSPAN  8-ft.  tran¬ 
sonic  wind  tunnel.  The  grid  and  carriage  data 
were  6%,  while  the  freestream  data  at  both  6% 
and  22%  were  available.  As  was  seen  in  Figures 
5-8,  the  wind  tunnel  data  were  in  excellent 
agreement  with  the  flight  test  results. 

CFD  Research  Corporation. 

L.  Hall11  presented  the  results  of  the  CFD  re¬ 
search  code  for  the  F-18/JDAM  configuration. 
These  results  were  significantly  different  from 
the  other  seven  codes  presented,  since  the  trajec¬ 
tory  calculations  were  run  in  a  time  dependent 
mode.  At  the  time  of  the  meeting,  the  trajectory 
had  run  for  only  .05  seconds;  however,  the  pre¬ 
dictions  shown  were  in  good  agreement  with  the 
test  data.  One  drawback  of  using  time  dependent 
(as  opposed  to  steady  state)  trajectory  calcula¬ 
tions  is  that  it  takes  a  very  long  time  to  get  one 
solution,  and,  if  any  of  the  parameters  change, 
another  solution  would  take  just  as  long. 

Air  Force  Wright  Research  Lab 

Tomaro  presented  the  F/A-18C/JDAM  carriage 
loads  and  trajectory  analyses  conducted  by 
AFRL/VAAC.  The  study  included  the  use  of 
computational  aerodynamic  (CFD)  and  6DOF 
rigid-body,  trajectory  generation  techniques.  The 
two  methods  were  not  coupled  into  a  single 
simulation  package:  the  CFD  method  simply 
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provided  the  aerodynamic  loads  database  to  the 
trajectory  generator  which  was  run  independently 
in  a  carriage-loads  decay  manner.  The  CFD  por¬ 
tion  of  these  analyses  used  the  AFRL/VAAC 
'Cobalt'  flow  solver. 

For  the  F/A-18C/JDAM  carriage  cases,  sev¬ 
eral  tetrahedral  grids  were  produced  with  the 
NASA  GridTool/VGRIDns  unstructured  mesh 
generation  system.  The  primary  grid  used  in  sub¬ 
sequent  carriage  loads  analyses  consisted  of  6.62 
million  cells  (half-model  for  symmetry)  with  vis¬ 
cous  boundary  layers  (approx.  4  million  cells) 
about  all  components  of  the  F/A-18C  and  JDAM. 
About  one-month  of  calendar  time  was  required 
to  generate  an  'appropriate'  mesh,  i.e.  no  negative 
volumes,  crossed  faces,  etc.  Subsequent  Cobalt 
solutions  required  about  10.1  GBytes  of  main 
memory  and  the  following  timings: 

M  =  0.962  M=  1.055 

50  nodes  (CPU's)  IBM  SP2  32  nodes  IBM  SP2 
17.69  hours  wall-clock  26.87  hours  wall-clock 
17.22  hours  CPU/node  6.27  hours  CPU/node 
(861  total  CPU  hours)  (841  total  CPU  hours) 

Isolated  freestream  JDAM  viscous  grid  gen¬ 
eration  and  flow  solutions  (alpha  and  beta 
sweeps)  required  less  than  3  weeks  turnaround 
for  both  Mach  numbers.  Thus,  within  two 
months  calendar  time,  the  CFD  portion  of  the 
Challenge  was  completed. 

After  the  carriage  and  isolated,  freestream 
JDAM  aerodynamics  were  provided  by  Cobalt, 
trajectories  were  generated  using  the  NAWCAD 
NAVSEP2  program.  The  carriage-loads  decay 
method  was  used  to  account  for  mutual  interfer¬ 
ence  effects  between  the  F/A-18C  and  JDAM.  In 
addition  to  aerodynamic  forces  and  moments, 
NAVSEP  requires  the  JDAM  inertial  properties, 
damping  coefficients,  ejector-model  characteris¬ 
tics,  and  then  it  calculates  trajectories  in  a  matter 
of  seconds  on  any  computer  platform. 

The  predictions  for  the  pitch  and  yaw  mo¬ 
tion  of  the  store  for  the  M  =  0.962  and  M  =  1.05 
were  in  excellent  agreement  with  the  test  data, 
although  the  pitch  attitude  was  somewhat  over¬ 
predicted.  This  implies  that  the  predicted 


carriage  pitching  moment  was  larger  than  that  in 
flight.  The  roll  motion  was  not  well  predicted. 
However,  since  rolling  motion  has  traditionally 
been  the  hardest  part  of  the  trajectory  to  predict, 
and  generally  has  little  influence  on  store  miss 
distances,  the  lack  of  rolling  motion  correlation  is 
of  small  consequence. 

Naval  Air  Warfare  Center. 

Woodson19  presented  comparisons  for 
SPLITFLOW,  USM3D,  and  PUMA,  an  unstruc¬ 
tured,  viscous  code  developed  at  NAWCAD. 

Only  viscous-store  results  were  computed 
using  the  USM3D  code.  The  code  for  both  cases 
was  run  for  2000  iterations  using  a  CFL  of  0.1 
initially  ramping  up  to  100  over  the  first  500  it¬ 
erations  and  then  continued  at  100  for  the 
remaining  1500  iterations.  The  solutions  con¬ 
verged  to  a  steady  state  value  in  approximately 
500  iterations  with  the  residual  reduced  about 
three  orders  of  magnitude.  The  solutions  were 
run  on  a  Cray  C90  and  required  315  MW  of 
memory  and  a  total  of  48.44  hours  of  CPU  time 
for  case  1  and  57.46  hours  for  case  2.  Multitask¬ 
ing  was  employed  using  ten  processors  for  a  wall 
clock  time  of  approximately  six  hours  for  case  1 
and  eight  hours  for  case  2  (average  concurrent 
CPUs  =  7.5  and  6.96,  respectively). 

Two  solution  approaches  using 
SPLITFLOW  were  conducted:  (1)  inviscid,  and 
(2)  viscous  around  the  JDAM  store.  By  assigning 
different  material  numbers  to  the  various  compo¬ 
nents  of  the  configuration  (i.e.  wings,  pylons, 
stores,  etc.)  different  boundary  conditions  may  be 
applied  so  that  the  prismatic  grids  may  be  gener¬ 
ated  only  on  those  parts  of  the  geometry  where 
viscous  effects  are  anticipated  to  be  important 
and  neglected  elsewhere.  Both  cases  were  run  for 
2000  iterations  using  a  CFL  number  of  1.0  and  a 
turbulent  CFL  number  of  0.1.  The  solutions  were 
run  on  a  Cray  C90  requiring  256  MW  of  memory 
and  a  total  of  58.56  CPU  hours  for  case  1  and 
81.29  hours  for  case  2.  Multitasking  was  em¬ 
ployed  using  four  processors  for  a  wall  clock 
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time  of  approximately  34  and  48  hours  (average 
concurrent  CPUs  =  1.75  and  1.7,  respectively). 
The  longer  run  time  for  case  2  was  caused  by 
sliver  cell  problems  aft  of  the  shock  at  the  pylon 
trailing  edge  so  a  smaller  global  cell  size  was  em¬ 
ployed  which  resulted  in  the  code  reaching  its 
maximum  number  of  cells  much  sooner  than  it 
did  for  case  1.  Both  cases  achieved  about  three 
orders  of  magnitude  reduction  of  the  residuals. 
The  inviscid  solutions  required  about  one  half  the 
run  times  of  the  viscous  store  results  (26.50  and 
40.56  CPU  hours,  respectively). 

Both  of  the  PUMA  runs  were  inviscid  and 
were  performed  using  between  32  and  64  nodes 
of  the  IBM  SP-2.  Each  run  was  converged  two 
orders  of  magnitude  using  first  order  spatial  ac¬ 
curacy.  Then  this  first  order  solution  was  used  to 
initialize  the  second  order  runs.  Unfortunately,  it 
was  possible  to  converge  the  second  order  solu¬ 
tion  only  about  1-2  orders  in  the  residuals. 

Although  the  results  presented  were  impres¬ 
sive  and  showed  good  correlation  with  the  test 
data,  the  JDAM  strakes  were  not  modeled  in  the 
analysis.  It  is  not  known  whether  the  results  ob¬ 
tained  were  fortuitous,  or  that  the  strakes  have 
little  impact  on  the  trajectories. 

Lockheed  Martin 

Welterlen20  presented  both  viscous  and  inviscid 
SPLITFLOW  results.  The  viscous  grid  required 
800,000  Cartesian  Cells,  and  an  additional 
1,044,207  prismatic  cells  on  the  surface.  The  in¬ 
viscid  solution  required  approximately  120  CPU 
hours  on  a  HP  9000  using  8  V-2250  processors. 
The  viscous  results  required  about  250  CPU 
hours  on  a  Cray  J-90. 

The  predicted  carriage  loads  were  in  close 
agreement  with  the  test  data;  however,  the  invis¬ 
cid  results  matched  the  carriage  loads  better  than 
did  the  viscid.  Since  the  SPLITFLOW  code  was 
coupled  to  a  six-degree-of-freedom  code  that  was 
developed  for  this  purpose,  it  was  not  possible  to 
determine  whether  the  relatively  poor  predicted 


trajectories  were  due  to  SPLITFLOW  or  the  tra¬ 
jectory  code. 

Aerosoft  Inc 

The  sixth  paper  from  Aerosoft  Inc  using  the 
GUST  solver  package  was  withdrawn.  McGroy 
later  presented  results  that  were  similar  to  the 
others  shown  in  this  paper,  but  these  have  not 
been  made  available  at  this  time. 

DSTO  Australia 

Fairlie  then  presented  computations22  using  the 
RAMPANT  code  supplied  by  Fluent  Inc. 

The  trajectories  followed  by  the  JDAM  after 
its  release  were  simulated  using  the  Defence  Sci¬ 
ence  and  Technology  Organization  Release 
Evaluation  Suite  (DSTORES).  This  approach  is 
similar  to  that  used  by  Tomaro  and  Woodson. 

The  initial  grid  as  input  to  the  RAMPANT 
Solver  consisted  of  just  over  1.05  million  tetra- 
hedra.  A  typical  RAMPANT  run  consisted  of 
about  200  iterations  with  the  value  of  the  CFL 
number  set  to  0.5,  followed  by  about  300  itera¬ 
tions  with  CFL  set  to  1.0.  This  was  generally 
sufficient  to  reduce  the  normalized  residuals  of 
continuity,  x-,  y-  and  z-momentum  and  energy  by 
between  two  and  one  half  and  three  orders  of 
magnitude  compared  with  their  initial  values.  At 
this  stage,  the  grid  was  adapted  in  regions  in 
which  the  static  pressure  gradient  exceeded  a 
particular  value  (initially  set  to  10\%  of  its 
maximum  value,  but  varied  depending  on  the 
number  of  tetrahedra  generated  in  the  new  grid) 
in  an  attempt  to  better  define  shock  waves.  After 
adaptation,  the  grid  generally  contained  some¬ 
what  more  than  1.25  million  tetrahedra.  The 
solution  was  then  iterated  for  up  to  a  further  500 
iterations.  After  an  initial  transient  created  by  the 
adaptation  of  the  grid,  the  residuals  once  again 
fell  to  their  previous,  or  even  lower  values. 

All  computations  were  carried  out  on  a  Sili¬ 
con  Graphics  Origin  2000  server.  This  machine 
has  sixteen  R 10000  processors  running  at 
250~MHz  and  is  equipped  with  4~GB  of  mem- 
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ory.  Typically,  the  initial  500  iterations  on  the 
un-adapted  grid  required  a  little  less  than  40 
hours  of  CPU  time  on  a  single  processor,  and  oc¬ 
cupied  approximately  460-MB  of  memory.  After 
adaptation,  the  additional  500  iterations  used 
somewhat  more  resources,  the  exact  amount  de¬ 
pending  on  the  number  of  tetrahedra  in  the 
adapted  grid.  While  the  vast  majority  of  the  cal¬ 
culations  presented  were  carried  out  on  a  single 
CPU,  the  RAMPANT  code  may  be  run  in  paral¬ 
lel.  Thus  far,  no  more  than  four  parallel 
processors  have  been  used,  yielding  a  speed-up 
of  just  over  3.8  compared  with  a  single  proces¬ 
sor. 

The  RAMPANT  predicted  carriage  loads 
were  input  into  the  Australian  six-degree-of- 
freedom  code,  in  conjunction  with  the  experi¬ 
mental  store  freestream  data,  and  the  trajectories 
were  calculated  in  a  manner  similar  to  those  in 
Reference  17.  The  Australian  code  has  an  un¬ 
usual  feature  that  allows  it  to  calculate  a  yaw 
restraint  between  the  pistons  and  the  store  during 
the  ejector  stroke,  which  lasted  for  approximately 
0.07  seconds.  The  yaw  attitude  would  have  been 
in  excellent  agreement  with  the  flight  test  data  if 
the  prediction  were  displaced  by  0.07  sec.  It  ap¬ 
pears  that  for  this  case,  the  trajectory  code  did  not 
properly  account  for  the  constraint  between  the 
pistons  and  store.  However,  there  have  been 
numerous  flight  test  cases  where  the  store  was 
clearly  constrained  in  yaw  during  the  ejection 
stroke.  The  constraint  feature  will  become  more 
useful  once  it  has  been  calibrated  with  flight  re¬ 
sults. 

The  pitch  attitude  prediction  was  in  reason¬ 
able  agreement  with  the  test  data,  although  it 
overpredicted  the  test  data  by  approximately  20% 
at  both  Mach  numbers.  The  roll  attitude  was  in 
excellent  agreement  with  the  flight  test  results; 
the  yaw  constraint  might  have  fortuitously  helped 
to  constrain  the  roll. 


NRCIAR  Canada 

The  last  paper  described23  the  quasi-steady  CFD 
approach  developed  at  the  Institute  for  Aerospace 
Research  (IAR)  of  the  National  Research  Council 
of  Canada  (NRC).  It  consists  of  three  different 
modules: 

1)  A  steady-state  3D  unstructured  inviscid  solver, 
FJ3SOLV. 

2)  A  6-DOF  Store  Separation  Model  (SSM),  and 

3)  A  grid  motion  technique. 

Each  of  these  modules  could  be  used  sepa¬ 
rately. 

To  apply  the  IAR  approach  to  the  F/A-18C 
JDAM  CFD  Challenge,  the  three  modules  were 
coupled  in  a  quasi-steady  mode  using  the  fol¬ 
lowing  methodology: 

For  a  given  store  position,  compute  the 
steady-state  aerodynamic  loads  acting  on  the 
store  using  FJ3SOLV.  Feed  the  CFD  predicted 
aerodynamic  loads  into  the  6-DOF  SSM  and,  for 
a  small  time  increment  compute  the  new  store 
CG  location  and  angular  orientations. 

If  grid  motion  is  possible,  move  the  store 
and  grid  nodes  using  spring  analogy  technique  to 
store’s  new  position  and  go  to  step  1.  If  grid  mo¬ 
tion  is  not  possible  or  grid  cells  become 
inadequate  after  node  movement,  move  the  store 
to  its  new  position,  generate  a  new  grid,  interpo¬ 
late  the  solution  to  the  new  grid  from  the 
previous  one  and  go  to  step  1.  Starting  from  a 
clean  configuration,  suitable  for  girding,  and  with 
a  time  step  of  0.02  sec.,  it  took  about  two  weeks 
to  compute  the  JDAM  trajectory  for  a  time  period 
of  0.24  seconds.  The  Canadian  results  were  very 
similar  to  the  others  shown. 

OTHER  RESULTS 

Two  other  organizations  that  tried  to  take  part  in 
the  Challenge  were  not  able  to  present  their  re¬ 
sults  in  Reno.  Their  approach  differed  from  all 
the  previous  papers  since  they  used  a  structured 
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grid  approach  based  on  the  chimera  formulation. 
The  results  presented  are  described  below. 

AFSEO 

The  Applied  Computation  Fluid  Dynamics 
(ACFD)  group  within  the  Air  Force  Seek  Eagle 
Office  (AFSEO)  also  computed  the 
F/A18C/JDAM  Challenge  cases24  using  a  fully 
time  accurate  CFD  simulation. 

The  ACFD  group  utilized  the  Beggar  flow 
solver  originally  developed  at  the  Air  Force 
Wright  Laboratory  at  Eglin  AFB  with  develop¬ 
ment  continuing  within  the  AFSEO  by  the  ACFD 
group. 

The  separation  of  the  JDAM  from  the 
F/A18C  at  Mach  0.962  was  simulated  using  the 
Beggar  code  assuming  inviscid  flow.  The  grid 
system  for  the  F/A18C,  the  JDAM,  and  assorted 
auxiliary  grids  utilized  a  total  of  39  single  block 
grids  and  12  multi-block  grids  with  a  total  of  95 
grids  and  2.8  million  grid  points.  The  JDAM 
grid  alone  contained  360,000  points.  The  grid 
system  was  generated  in  approximately  one  man 
month. 

The  separation  simulation  was  run  at  a 
physical  time  step  of  one  millisecond  and  was 
terminated  at  a  solution  physical  time  of  0.42 
seconds.  The  solution  was  run  on  16  processors 
for  the  flow  field  solution  with  another  2  proces¬ 
sors  used  for  the  grid  assembly  process.  The 
execution  time  varied  somewhat  with  an  average 
wall  clock  time  of  160  seconds  per  time  step  on 
an  SGI  Origin  2000  with  250MHZ  MIPS  R 10000 
processors.  Thus,  a  simulation  out  to  0.25  sec¬ 
ond  in  physical  time  could  be  performed  in  less 
than  12  hours.  Each  flow  solver  process  utilized 
between  75  and  132  Megabytes  of  memory  with 
the  total  memory  requirements  of  1.6GB  for  the 
flow  solver.  Each  of  the  two  grid  assembly  proc¬ 
esses  utilized  another  512MB  of  memory. 

The  agreement  shown  in  general  was  excel¬ 
lent,  with  the  inviscid  results  slightly  over 
predicting  the  pitch  and  yaw  angles.  The  prelimi¬ 
nary  viscous  results  showed  an  improvement  in 


the  agreement  with  the  flight  test  data.  The  invis¬ 
cid  prediction  of  the  roll  angle  was  generally 
good  and  captured  the  general  trends.  The  vis¬ 
cous  agreement  degraded  at  later  times  when  the 
JDAM  was  at  large  yaw  and  pitch  angles. 

AEDC 

AEDC  performed  time-accurate  viscous  compu¬ 
tations  to  simulate  the  trajectory  simulations  for 
both  flight  release  conditions25  utilizing  the  chi¬ 
mera  overset  grid  approach.  The  process  for 
predicting  time-accurate  body  motion  relies  on 
four  codes,  NXAIR  to  solve  the  fluid  dynamic 
equations,  PEGSUS  to  define  the  inter-grid 
communications,  FOMOCO  to  compute  the  store 
loads  for  overlapping  surfaces  entities,  and 
SIXDOF  to  solve  the  rigid-body  equations  of 
motion. 

The  F-18C  surface  definition  was  prepared 
from  a  CAD  definition.  A  significant  portion  of 
the  effort  involved  preparing  the  aircraft  surface 
from  the  CAD  definition,  which  was  accom¬ 
plished  in  approximately  one  month.  Volume 
grid  generation  and  setting  up  the  PEGSUS  in¬ 
puts  required  about  two  to  three  weeks.  To 
reduce  the  number  of  grid  points  to  define  the 
boundary  layer,  wall  functions  were  utilized.  The 
overall  grid  system,  comprised  of  7.0  x  106  mesh 
points,  is  distributed  over  66  individual  overset 
meshes  (5.2  x  106  points  over  47  meshes  for  the 
F-18C  and  1.8  x  106  over  23  meshes  for  the 
JDAM).  All  detail  of  the  JDAM  was  modeled 
including  the  strakes  and  fin  gaps.  The  aircraft 
engine  duct  was  modeled  to  compute  flow 
through  the  duct.  Lateral  symmetry  about  the 
aircraft  center  plane  was  assumed,  and  only  the 
port  side  of  the  aircraft  was  modeled. 

Given  the  release  conditions  and  ejector 
model,  the  two  separation  trajectories  were 
simulated  by  using  the  aforementioned  codes. 
The  turbulent  Navier-Stokes  equations  were 
solved  with  the  two-equation  SST  turbulence 
model.  Duct  flow  was  established  to  a  corrected 
mass  flow  rate  at  approximately  145  lbm/sec  for 
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both  release  conditions.  The  steady-state  solu¬ 
tions  of  the  flow  field  about  the  carriage 
configuration  were  performed  until  convergence 
was  achieved  for  the  store  loads  to  approximately 
three  decimal  places.  The  steady  state  and  tra¬ 
jectory  computations  were  performed  on  a  SGI 
Origin  2000  R 10000  and  level  loaded  over  16 
processors.  Computations  to  determine  the 
steady-state  carriage  loads  required  approxi¬ 
mately  600  steps  for  each  case.  The  time-accurate 
computations  took  500  time  steps  to  compute  the 
0.4-second  trajectory.  The  total  CPU  time  to 
complete  one  case  (including  both  the  steady 
state  and  dynamic  portion  of  the  computations) 
was  2900  CPU  hours  (1400  for  the  steady-state 
solution  and  1500  for  the  dynamic  solution). 
Approximately  25  percent  of  the  time  in  the 
time-accurate  portion  of  the  problem  are  required 
by  PEGSUS.  Because  computer  resources  had 
to  be  shared  with  other  users,  only  part-time  us¬ 
age  of  16  processors  was  available  and  the  wall 
clock  time  to  complete  each  trajectory  simulation 
(steady  state  and  dynamic  portion)  was  two 
weeks.  With  dedicated  usage  of  64  processors  on 
an  Origin  2000,  the  computations  could  be  com¬ 
pleted  in  less  than  2  days. 

Comparisons  between  computed  orientation 
and  flight  telemetry  data  showed  excellent 
agreement  in  pitch  and  roll  while  the  computed 
yaw  showed  a  slightly  larger  nose  outboard  an¬ 
gle.  The  computed  store  displacements  and  miss 
distance  for  the  supersonic  case  showed  some¬ 
what  better  agreement  than  for  the  subsonic  case. 

Summary  of  ACFD  Challenge  II 

The  quality  of  the  invited  papers  and  presenta¬ 
tions  reinforced  the  approach  used  by  the  AFCD 
Challenge  sponsors.  However,  taking  these 
presentations  as  representative  of  state  of  the  art 
for  applying  current  CFD-based  tools  for  stores 
carriage  and  separations  indicates  that  wind  tun¬ 
nels  will  still  be  relied  on  for  the  provision  of  the 
major  part  of  the  aerodynamic  data  on  which 


stores  certification  are  to  be  safely  based.  Indeed 
it  is  acknowledged  that  the  CFD  solutions  were 
in  the  majority  of  cases  within  the  error  range  of 
the  wind  tunnel  and  flight  test  data.  Accuracy 
would  not  therefore  seem  to  be  issue,  but  rather 
the  time  required  to  produce  a  solution  needs  to 
be  decreased  significantly.  Given  this  develop¬ 
ment  CFD-based  tools  should  become  far  more 
prevalent  in  use  during  Requirements  Definition 
and  Systems  Engineering  trade-off  studies  for  the 
aircraft  and  stores  thereby  reducing  the  likely 
hood  of  expensive  aircraft  and/or  store  redesign 
after  hardware  has  been  made. 

One  other  general  result  was  the  consensus 
that  improvements  in  the  ejector  modeling  and 
ejector  foot/store  interaction  during  the  ejection 
needed  to  be  accomplished. 

One  of  the  principal  drawback  of  CFD 
Challenge  II  was  that  all  the  CFD  results,  using 
both  Euler  and  Navier  Stokes,  as  well  as  a  simu¬ 
lation  that  ignored  the  JDAM  canards  gave 
similar  results.  Does  that  mean  that  Navier 
Stokes  formulation  does  not  have  to  be  used,  or 
were  the  test  cases  selected  fortuitous  for  the  in- 
viscid  formulation.  Indeed,  Welterlen  showed 
that  his  inviscid  calculation  was  superior  to  the 
viscous  one.  Since  diagnostic  data  were  not 
available,  it  is  impossible  to  say  whether  the 
SPLITFLOW  viscous  formulation  was  at  fault,  or 
that  the  inviscid  results  had  a  fortuitous  canceling 
error.  It  was  the  consensus  of  the  participants 
that  another  CFD  Challenge,  one  that  would  have 
diagnostic  data  (store  and  wing  pressures)  was 
merited. 

FUTURE  PLANS 

Recently,  It  has  been  shown  that  (pressure  sensi¬ 
tive  paint)  PSP  can  be  used  to  evaluate  store 
pressures;  furthermore,  wind  tunnel  test  data  on 
the  F-18  aircraft  with  and  without  the  TFLIR 
(targeting  forward  looking  infrared  radar)  have 
demonstrated  extreme  Mach  number  sensitivity. 
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as  may  be  seen  in  Figure  9.  This  effect  was  con¬ 
firmed  by  recent  F-18/MK-84  flight  test  results. 

It  is  anticipated  that  ACFD  Challenge  III 
will  make  use  of  PSP  wind  tunnel  results  for  the 
F-18  aircraft  with  the  TFLIR  on,  as  well  as  the  F- 
18/MK-84  flight  test  data. 

JDAM  MK-83  STATION  3 


FIGURE  9  Mach  Sweep 
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